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Abstract. A generalised stochastic Gross-Pitaevskii equation describing a partially 
condensed trapped Bose gas with rotating thermal component is presented. We elucidate 
the manner in which the rotation changes the role of the high energy cutoff and introduces 
centrifugal effects in the classical field evolution. The rotation of the cloud means that Bose- 
enhanced collision processes occur preferentially into states with the same angular momentum 
as the thermal cloud, thus favouring vortex formation. We use the formalism to obtain a first 
principles theory of vortex lattice formation caused by thermodynamic instability. 
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1. Introduction 

1.1. Background 

The experimental achievement of dilute gaseous Bose-Einstein condensation [1] began a new 
era of intense interest in atomic, molecular and optical physics. Among their many fascinating 
properties, trapped dilute gas Bose-Einstein condensates (BECs) exhibit the properties of 
superfluidity in a precisely controllable and theoretically quantifiable manner [2]. At very 
low temperatures the trapped gas behaves like an ideal superfluid described by a macroscopic 
order parameter. Two persistent themes throughout the field are the Gross-Pitaevskii 
(GP) description of the condensate [3-5], and the Bogoliubov theory of its quasiparticle 
excitations [6]. The GP theory is an approximate zero temperature description of atomic BEC 
that accurately predicts a wide range of BEC behaviour, including ground states, dynamics, 
vortices and solitons [2]. However, it neglects some important quantum mechanical effects 
which we discuss below. The first quantum correction to the classical field theory of Gross 
and Pitaevskiif is provided by the Bogoliubov theory which gives a quantum description of 
weak excitations near zero temperature - the quasiparticle excitations of the condensate. 

1.2. Finite temperature effects and vortex formation mechanisms 

However, there are interesting dynamics that require a high temperature theory for their 
description, which necessarily must go beyond the GP and Bogoliubov theories. By high 
temperature we mean that the temperature is significant compared to the single particle 
energies of the trapping potential. For a harmonic potential with angular frequency lu the 
condition is 

hLu^ksT, (1) 

where fc^ is Boltzmann's constant, implying that thermal fluctuations will play an important 
role in the system. The process of vortex formation in BEC involves aspects of thermal, 
nonlinear and many body quantum physics, and our aim is to provide a theory that correctly 
represents the important physics of each of these. 

1.2.1. Stirring a condensate Several experimental groups have generated vortices in a 
trapped BEC using a rotating elliptical perturbing potential effected with a focused detuned 
laser beam [7-13]. The physical mechanism is the introduction of a symmetry breaking 
perturber which couples the ground and excited states of the condensate. Perturbative analyses 
of stirring [14-17] have found the critical rotational speed at which dynamical instability sets 
in leading to growth of the vortex state, and shown that the critical frequency is given by the 
angular form of the Landau critical velocity for superfluid dissipation. Thus it would appear 
that dissipation is essential to the process, as it is required for the system to surmount the 
energy barrier which otherwise prevents vortex penetration from the condensate edge [18]. 

In a non-perturbative simulation of a localized stirrer based on the Gross-Pitaevskii 
equation [19], it was found that a kind of nonlinear Rabi cycling occurs where the vortex 
periodically cycles from infinity into the condensate interior In the absence of dissipation a 
stable lattice does not form. The physical role of the stirrer in experiments, therefore, must be 
both to excite multi-pole modes of the condensate, and to create a rotating thermal component 
which can provide the requisite dissipation. 

f The use of the word 'classical' here requires some justification; a full discussion is given in section l43l 
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The physical mechanism arises from the fact that a condensate containing a vortex is 
thermodynamically unstable in the laboratory frame, so that the system energy can be lowered 
by irreversible excitation of the anomalous mode [20] . In the presence of dissipation the vortex 
spirals out of the condensate. However, for this to occur the dissipation must arise in the 
lab frame, corresponding physically to a condensate immersed in a stationary thermal cloud. 
Conversely, if the cloud is rotating faster than the critical frequency for vortex nucleation, 
the lowest energy state in the co-rotating frame is a condensate containing a stable vortex 
lattice [18]. 

7.2.2. Cooling a rotating cloud It has long been appreciated that phase defects can 
be trapped during a phase transition [21-23], thus providing another vortex formation 
mechanism. This has been demonstrated theoretically for Bose-Einstein condensation by 
Anglin and Zurek [24], Marshall et al [25] and Davis et al [26], but the mechanism has not 
yet been observed in BEC experiments. 

In a series of elegant experiments, Haljan et al [27] used a novel technique to create 
large vortex lattices. Unlike the other nucleation experiments at that time which acted on 
a pre-formed condensate with some kind of mechanical stirrer, this experiment involved 
evaporatively cooling a rotating thermal cloud to quantum degeneracy. In this way, it became 
possible to directly nucleate a condensate containing a disordered array of vortices with the 
same sense of circulation as the thermal cloud; nevertheless, direct evidence that this process 
involves defect pinning during the phase transition has not been forthcoming. 

Thus far there has been no theory of the formation of vortex lattices that includes the 
effects of thermal and quantum fluctuations, and, perhaps more importantly, is capable of 
describing the rotating condensation experiments of Haljan et al [27]. The development of 
such a theory has, in large part, been the motivation for this work. 

1.2.3. Including dissipation One theoretical approach to studying formation dynamics has 
been to introduce some form of phenomenological damping term in the GPE [28,29], which 
is intended to model the effect of a rotating thermal cloud. In the absence of an exact 
microscopic derivation several different approaches have been taken as to the form such a 
damping term should take. Tsubota et al [29] introduced a small imaginary coefficient in 
the time derivative of the GPE, in the rotating frame of the thermal cloud. Gardiner et al 
[30] derived a similar equation, but from a more physical point of view, by replacing the 
condensate chemical potential that arises in the GPE from a local energy conserving kinetic 
theory with a time derivative, leading to the Gardiner- Anglin-Fudge equation. Both of these 
approaches lead to the same physics as the form of the vortex growth equation we derive in 
section0of this paperf , which is found from a clear microscopic derivation within controlled 
approximations. 

7.2.4. Hamiltonian chaos, Kolmogorov turbulence and the Wigner representation A paradox 
arises when considering the mechanism of vortex formation via mechanical stirring. On one 
hand, it seems clear that some kind of dissipation is required which would be provided by 
a rotating thermal component. On the other hand, time evolution in quantum mechanics is 
unitary, so that if we were able to fully simulate the quantum many body problem for stirring 
a BEC initially in a ground state of the trap in the laboratory frame, we would expect to find 
that the system never reaches equilibrium in the rotating frame of the stirrer. 

f The exact equivalence arises in tlie limit of low damping. 
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Putting aside the origin of dissipation in quantum mechanics which can often be trivially 
resolved by recourse to the law of large numbers, a similar paradox arises in the Gross- 
Pitaevskii theory, which is also a Hamiltonian theory. In practice one observes that the 
dynamics are irreversible - non-equiUbrium initial states will evolve into quasi- stationary 
states with ergodic properties, wherein time evolution generates a thermal ensemble [31]. An 
important feature of the GPE is that the motion can exhibit features of classical Hamiltonian 
chaos [32]. It is therefore possible that a kind of internal dissipation can provide a mechanism 
for vortex nucleation under the right circumstances. This possibility has been considered in 
the work of Lobo et al [33] who used the classical field method to study the effect of a rotating 
drive on a BEC. The authors found that in the absence of any external dissipation or initial 
noise the nonlinear turbulence of the classical field evolution in the vicinity of a dynamical 
instability [34] was sufficient to relax the classical field toward equilibrium - although the 
relaxation occurs over a rather long timescale. 

The resolution of this paradox comes from the fact that the classical field evolution of 
the GPE does not provide a complete physical description of the many body quantum system 
of interest. In fact, we can expect the Kohnogorov spectrum of turbulence exhibited by the 
GPE near the instability [35] to be associated with the generation of a physical component 
of the gas that cannot be described within GP theory - a thermal component consisting of 
many weakly occupied modes. In Kolmogorov turbulence energy cascades from long to 
short wavelengths, and it is the turbulent dynamics at short wavelengths which are difficult 
to observe, difficult to simulate accurately within GP theory, and the reason why the GP 
evolution leads to dissipation. 

Indeed, the method of dynamical evolution of the GPE with initial noise arises from 
a Wigner representation of the quantum field, provided that terms beyond the GPE in the 
evolution which account for important quantum mechanical effects are truncated [36]. The 
truncation is usually justified if all of the modes included in the GP description are highly 
occupied], which is clearly not the case during a dynamical instability since many high 
energy modes are excited. Thus, the physical resolution of the paradox must come about 
from a theory that accurately describes the thermal cloud, and in particular its formation and 
influence on the coherent component of the gas. 

There is a distinct need, then, for a theory that accurately describes the rotating thermal 
component of a trapped Bose gas, and which includes the stimulated collisions between the 
condensate and noncondensate atoms. In this work we do not address the formation of the 
thermal cloud, but we instead go some way in the direction of a theory of vortex formation 
by providing a theory of the influence of the thermal cloud, a description which is necessarily 
stochastic in its formulation. 

1.3. SGPE theory 

One particularly successful way of handling both finite temperature and many body effects 
consists of treating the trapped Bose gas as an open system. The finite temperature Bose gas 
treated as a condensate coupled to a thermal reservoir which generates dissipative effects. The 
first derivation of a generalised stochastic Gross-Pitaevskii equation (SGPE) from this point 
of view for trapped BEC was provided by Stoof [39, 40] using functional techniques, starting 
from a Keldysh path integral formulation of Quantum Field Theory. The differences between 
the theory of [39] and the present SGPE theory are discussed in [41], and we refer the reader 
to this paper for detailed discussion; however, the most notable difference is the absence of a 

t Recent work by Nome et al [37] shows that it is possible to relax the high mode occupation criteria somewhat, in 
particular, for the condensate collision scenario of [38] where many vacuum modes are included in the description. 
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particular kind of interaction between the condensate and thermal cloud, which we refer to as 
the scattering term, which will be described in section l33l This term does not directly change 
the occupation of the condensate band, and so often its effect can be justifiably neglected. 
It does however couple to superfluid flows in the condensate band, so that, for example, 
we may expect it to play an important role in vortex motion. Stoof 's approach has certain 
advantages, most notably, since it is based on a Lagrangian formulation it is amenable to 
variational techniques, and a Gaussian variational method has been developed to great effect 
by Duine et al [42,43]. Nevertheless, as we have noted, there are certain effects which are 
difficult to capture using path integral methods, and for this reason we use the open systems 
approach developed in quantum optics [44,45]. 

In this paper we derive a generalised SGPE describing a rotating partially condensed 
trapped Bose gas, by adapting the theory developed by Gardiner et al [30,41]. In order to 
handle both a coherent condensed fraction and a high energy incoherent fraction, the SGPE 
formalism separates the partially condensed system into a low energy subspace of modes - 
the condensate band, and its high energy counterpart - the noncondensate band. The latter is 
treated as being thermalised, thus playing the role of a reservoir which damps the dynamics of 
the condensate band - in the rotating frame of the thermal cloud. An important feature of the 
approach, which also distinguishes it from that of [39], is that the condensate band contains 
both the condensate and its excitations. 

Using this approach we obtain stochastic equations of motion similar to those describing 
the nonrotating scenario of [41]. In the final analysis we show that the rotation generates 
additional centrifugal terms in the description of the condensate band and the thermal cloud. 
Since there are also some corrections to the derivation and final form supplied in [41], our 
approach in the present paper is to provide a reasonably self contained generalisation of that 
work and to supply a number of necessary corrections. 

1.4. Qualitative effect of the rotation 

There are two scenarios of interest in which a non-condensed Bose gas in some state of motion 
may come into thermal equilibrium. 

• In the first, the gas can have a constant center of mass motion in the laboratory, either 
confined to an atomic waveguide [46], or in a translating harmonic trap. It is well known 
that within classical statistics the thermodynamic properties of a uniformly translating 
gas are unaltered by the motion [47] - a consequence of Galilean invariance. This 
invariance trivially carries over to quantum statistics, and the dynamics of Bose-Einstein 
condensation for a uniformly translating gas are thus completely unaltered by the center 
of mass motion; the equations of motion are identical to those for the stationary system, 
modulo a shift in center of mass momentum, after a Galilean transformation into the 
moving frame. 

• The second situation of interest occurs when the gas is confined in a cylindrically 
symmetric trap, and has non-zero angular momentum. The gas is thus at rest in a frame 
rotating about the symmetry axis of the trap. In this case the motion of the cloud has a 
profound effect on the dynamics of Bose-Einstein condensation. 

A central feature of what follows is that the theory is developed in the frame co-rotating with 
the thermal cloud. Within classical statistics the only effect of rotation is to alter the dynamics 
of the thermal cloud such that it evolves in a centrifugally modified effective potential [47]. 
Within quantum statistics, however, there are important centrifugal terms which appear in the 
growth and scattering terms of the SGPE which we derive. It will be shown that the full theory 
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Figure 1. Schematic of the separation into condensate and noncondensate bands. States 
beneath the energy En, defined in the rotating frame, form the description of the condensate 
band which warrant a quantum mechanical description, and the states above Er form the 
noncondensate hand which consists of high energy, approximately thermal states. The S-wave 
scattering description of the collisions means that states with energies above fi^q^ /2m in the 
rotating frame are also eliminated from the theory [41]. 



of [41] can be used with the appropriate changes to the single particle Hamiltonian and the 
thermal distributions in the rotating frame. 

2. The system 

2.1. The cold-collision Hamiltonian 

In the cold-collision regime described by S-wave scattering length a [2] the second-quantised 
many body Hamiltonian for a dilute Bose gas confined by a trapping potential l^(x), 
formulated in the rotating frame defined by angular frequency Jl is [48] 

H^H,p + Hi (2) 

- where the single-particle Hamiltonian is 

H^p= / d'''x7/>t(x) f — — + ihn- (xxV) + V{x)j (3) 

and the interaction Hamiltonian is 

Hi = \j d'x ^t(x)V;t(x)V'(x)V'(x), (4) 

and where u = Anfi^a/ni, m is the atomic mass and a is the S-wave scattering length. 
As shown schematically in figure [0 the states of the trapped system are divided into 
the condensate band of states with energy beneath the cut-off Ep, and the remaining 
noncondensate band of high energy states. The separation into distinct bands allows the 
two regions to be treated using different techniques. The high energy states are to be treated 
using the master equation techniques of quantum optics, and will physically play the role 
of a thermal reservoir for the low energy condensate band. The condensate band modes are 
assumed to be, at least, moderately occupied so that some of the quantum mechanical features 
of the condensate band may be consistently neglected. The field operator is 

^(x) = </)(x) + Vnc(x) = Vi^i^) + QV(x), (5) 
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where the noncondensate field i/'nc(x) = Q^(x) describes the high energy thermal modes, 
and the condensate band is described by the field (/'(x) = V^fj(x.). The the orthogonal 
projectors V + Q — 1 sae defined with respect to the single particle basis of the trap. The 
condensate band projector takes the form 

P = ^|n)(n| (6) 

n 

where the bar denotes the cut-off and n represents all quantum numbers required to index the 
eigenstates of the single particle Hamiltonian. In the spatial representation the eigenstates are 
denoted by y„(x) = (x|n) and the action of the projector on an arbitrary wavefunction x(x) 
is 



(7) 



2.2. Separation of the Hamiltonian 

Using (|5} the full Hamiltonian (|2j can be written as 

H = Hq + Hjc + Hnc (8) 

where Hq involves only condensate band operators, Hj^c involves only noncondensate band 
operators and 

Hi,c = H^i]l + H^i% + H^i% (9) 
denote interaction terms involving one, two or three condensate band operators: 

Ho = d^x (/.t (x) f — — + ihn • (x X V) + l/(x) j 

d3x0t(x)0t(x)0(x)0(x) (10) 



Hnc - / d^x Vjvc W ( — — + • (x X V) + y(x) ) VArc(x) 



2m 

d^X ?/']vc(x)'0]vc(^)V'WC'(x)V'Arc(x) (11) 

Hf). =ujd?x V]vc(x)V']vc(x)V'ivc(x)<^W + li.c. (12) 
Hf).=u f d3xV']vc(x)'/'ivc(x)./)t(x)0(x)+h.c. 



+ V'jvc WV-jvc W'/'(x)0(x) + h.c. (13) 

H^j^l = ujd^x <^t(x)0t(x)0(x)V'A,c(x) + h.c. (14) 

We note that there are additional terms arising from the single particle Hamiltonian involving 
one ipNc{x) and one 0(x) but these play no role in the following theory because the 
noncondensate band operators have no mean field. 
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3. Condensate band master equation 

3.1. Formal derivation of the master equation 

The master equation for the condensate band density operator 

PC = trjvc (p) (15) 
is found from the von Neumann equation for the total density operator 

- j^[Ha + Hi + Hncp] 

= {Cq + Ci + Cnc)p (16) 
by defining the density operator projector 

VcP = PNC '^i'CNc (p) ^ PNC PC, (17) 

together with its orthogonal projector Qc = 1 — Vc', one then uses standard methods [45] to 
find the equation of motion for the Laplace transform of v{t) = Vcp{t) 

sv{s)-v{Q) = {C^ + VcCi)v{s) 

+ VcCi[s ~Cq-Ci- Cnc]~^QcCiv{s). (18) 

At this point two approximations are made. Firstly the evolution is approximated by 
neglecting the interaction superoperator £/ in [ ]^^, which amounts to treating the interactions 
as a perturbation to the bare system evolution. After using the Laplace convolution theorem 
to invert ( I18> , a Markov approximation is made for the time integral by neglecting the time 
dependence of v{t) over the support of the integral and extending the upper limit of integration 
to infinity. The approximation is that the reservoir correlation time is small compared to the 
timescale of system evolution. The result is 

v{t) = {Co+VcCj)v{t) 

+ S^VcCiJ^ dTexp{{CQ + CNc)T}QcCi^v{t). (19) 
We now evaluate the terms in ( I19> . 

3.2. Hamiltonian terms 

The term {Cq + Vc^i) in ( I19> can be written in terms of an effective condensate band 
Hamiltonian 



He = Hq + iJforward (20) 



where 



-ff forward = 2u J UnC {^)4>H^)4'{^) (21) 

and the noncondensate density 

nArc(x) EE {ip\-x.)ip{:x.)) = tipfc {pNCipH^)^i'^)) (22) 
plays the role of an effective potential. The Hamiltonian term is 

^c|Ham = -|[^C,Pc]- (23) 
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Figure 2. Schematic of the processes arising from the interactions between the condensate 
and noncondensate bands. In (a) two noncondensate band atoms collide. The collision energy 
is transferred to one of the atoms with the other passing into the condensate band. In (b) a 
condensate band atom collides with a noncondensate band atom with no change in condensate 
band population. The time reverse of the two processes also occurs. Figure adapted from [49]. 



3.3. Interactions between condensate and noncondensate bands 

The master equation terms arising from evaluating (I19> for ijj^p and Hfjj are 

(^jvc(x')V'ivc(x')'<c(X:^)'^JVc(x,T))[,/.(x,T)pc,0^(x')] 

+ (V^j^c(x, r)Vjvc(x, T)^j^^(x')VArc(x')')[0^(x'), Pc0(x, r)] 

+ (^jvp(x') VjvclxOV-jvcC^' t)iAjvc(x, t)2) [0t (x, r)pc, 0(x')] 

+ (<c(x,r)V'ivc(x,T)2^]vp(x')Vjvc(x'))[0(x'),Pc0^(x,T)]} (24) 



(V-jvcC^' t)V'7Vc(x, T)V'jvc(^')V'JVc(x')) 

X [pc0t(x,r)0(x,r),0t(x')0(x')] 

+ (■*/'jvc(^')^Afc(x')V']vc(^' T)?AArc(x, t)) 

X [0t(x')0(x'),0t(x,T)^(x,T)pc]} (25) 

where the time evolved field operators are defined as 

0(x,t) ^e'"'^-*/'''4>{-K)e-~'"<^'/^- (26) 

V'jvclx, = e^^«^*/Vivc(x)e^*^"^*/^ (27) 

There are also terms that arise in \\9\ from a mixture of terms from H^pQ and Hj"^^, and from 

-ff) p alone. All of these contributions can be seen to be prohibited by energy and momentum 
conservation, which is to say they are non-resonant, and as such are very small. 

The physical processes corresponding to the Hj^^, (Growth) and Hj^^j (Scattering) terms 
are shown schematically in figure 121 The growth terms correspond to scattering between 
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condensate and noncondensate atoms whereby particle transfer can take place. The so- 
called scattering terms involve collisions between condensate and noncondensate atoms which 
nevertheless conserve the populations in each band. 

3.4. Reservoir correlation functions 

We now introduce the approximations for the correlation functions of the noncondensate band 
operators. 

3.4.1. Hartree-Fock factorisation For a thermal gas one may use Hartree-Fock factorization 
to expand high order moments in terms of second order moments. The correlation function 
MAI , for example, becomes 

(V']vc(x')V'a'c(x')V]vc(x,^)Vwc(x,t)) = 2(?A]Yc(x')?AArc(x,T)) 

X (V'jvc(x')V']vc(^'^)> 
X (V'wc(x')V']vc(^>^)) (28) 

where terms of the form (V']vp(x') V'jvc (x')) are neglected because they do not lead to energy 
conserving processes in the master equation. 

3.4.2. Semi-classical description of the thermal cloud in the rotating frame The single 
particle Wigner function for the noncondensate field is defined as 



^(x, K) = y d\ + v/2) V'JVC (x - v/2)) 



e*^". (29) 



In the semiclassical description, the noncondensate band Wigner function for a Bose gas 
confined in a cylindrically symmetric trap in equilibrium in the rotating frame defined by 
VL = Viz, with chemical potential /z, takes the form [47, 50] 

exp [(fia;(x, K) — ^)/kB-L \ — 1 



where 



nuj(x, K) = — hn-(xxK) + v(x) 

= ^{K-mnxx/hf + VM^), (31) 
Zm 



and we have introduced the effective potential 



V;.(x) . Vi.) - I^fci^ (32) 

^ m{ujl - yi^ fr^ _^ mujlz^ ^^^^ 



where r±_ = -\- y'^ is the distance from the symmetry axis of the trap. From J3H we see 
that, in terms of the rotating frame momentum, the rotating frame semi-classical distribution 
is a lab frame equilibrium distribution with a centrifugally modified effective potential. The 
Coriolis force thus has no effect on the equilibrium state of the rotating cloud within the 
semi-classical picture [47]. 
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Introducing the variables u = (x + x')/2, v = x' — x, we can write 

(V'jvc(x')V'ivc(x, r)) = (7Aj^p(u + v/2)Vjvc(u - v/2, r)) 

-TTTT^I d^Ki^(u,K)e-*^-^-''^("'^)^ (34) 

(^jvc(x')V']vc(^. ^)) « jJ-^ I '^'K [1 + F(u, K)]e^K.v+i<.(u,K). (35^ 

where we have made a semi-classical approximation for the time development of the 
noncondensate field operator and used the fact that in the semiclassical limit F{u, K) is only 
significant in the noncondensate region, which is given by the condition 

nj^{\i,K) > Er. (36) 

We note that the rotating frame cutoff condition is more complicated than the laboratory 
frame cutoff since the former is defined by a position dependent rotating frame momentum. 
However, the only dependence of F{u, K) on K — mJl x x/?i is implicit, through hio{u, K). 
In practice one can introduce an appropriate change to energy variables in the necessary 
integrals - effectively eliminating the position dependence, so that there is no difficulty in 
handling this complication. 

Equations ( 13 1> . ( I34> . and ( I35> form our description of the thermal cloud; the lab 
frame appearance of (|3T) means that much of the Quantum Kinetic Theory [51-53] that 
has been developed to describe condensate growth also applies in the rotating frame with 
straightforward modifications. However, due to the rotation there are some important changes 
to the condensate band description. 

3.5. Master equation 

3.5.1. Growth terms The terms in the master equation that lead to condensate growth are 
from the interaction Hamiltonian term H^P(j, corresponding to the master equation term \2A\ . 
This term will be rewritten using the approximate description of the reservoir correlation 
functions introduced above; further, we will neglect the principal value part of the time 
integrals in (|24j and (|25} 

1-0 

dr exp {—ilot) — 'k5{lj) + iP/oj w ■k5{lo), (37) 



which is physically consistent with the Markov approximation since it ensures conservation 
of energy during collisions between condensate and noncondensate atoms. 
The growth terms can be written in terms of the amplitudesf 

G<+'(u,v,e) = .^^-^ y^^Ki j Jd^Ks Sicj^ + uj2 ~ u;3 - e/h) 

(u, V, e) = j^^^ J c^'Ki J d^Ka J d^Kg S{lo, + u;2 ~ lus ~ e/h) 

X [l + i^(u,Ki)][l + F(u,K2)]F(u,K3)e-'(Ki+K^~K^)-", (39) 
and the condensate band operator 

Lcc^{-k) = [^{-x),Hc] (40) 

f This corrects an extra minus sign in tlie defining equation (56) of [41]. 
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Pc|g_,h= / / rf'v[{G<-'(u,v,Lc)0(u-v/2)}pc,0nu + v/2)] 

- j dPn J d^v [pc{G<-'(u,v,-Lc)</'Uu-v/2)},(^(u + v/2)] 
+ Jd^'u y"d3v[{G<+'(u,v,-Lc)0^(u-v/2)}pc,'/'(u + v/2)] 

- j d'nj d\ [pc {G<+' (u, V, Lc)<^(u - v/2)} ,cj>\n + v/2)] . (41) 

3.5.2. Scattering terms In terms of the scattering amplitude 

M(u, V, e) = j t^'Ki j d^Ka 5{uji - c^2 - e/^) 

X F(u,Ki)[l + F(u,K2)]e*(^i"^^)-"' (42) 

and the operator 

C/(x) =^t(x)^(x), (43) 
the scattering terms, given by \25\ . take the form| 

PcL.tt= j / rf'v[C/(u + v/2),pc{M(u,v,Lc)C/(u-v/2)}] 

+ J d^^J iUH^, V, -Lc)C/(u - v/2)} pc, C/(u + v/2)] (44) 

5.5.5. Forward-backward relations When the noncondensate band is described by the 
thermal equilibrium form given by (I30> . the growth and scattering amplitudes satisfy the 
forward-backward relations [41, 52] 

G'-'(u,v,e) = e("~'^)/'''«^G<+>(u,v,e), (45) 

and 

M(u, v, e) = e""/'=«^Af (u, v, -e). (46) 

3.5.4. Full master equation and its stationary solution Defining the condensate band 
number operator 

Nc = J d^Kcj)\x)(j)ix), (47) 
we can write the stationary solution of the full master equation 

PC = PC I Ham + Pc\g,,^tl. + PcL.tt ^^8) 

as 

^-"^"H )■ ^^'^ 

f This corrects a misprint in equation (59) of [41] wliere Lc appeared in place of —Lq in the second and third lines 
of 1411 . Consequently the identities 1501 - 1521 required to estabhsh the validity of the equilibrium solution appeared 
incorrectly in [41]. 

I Note the different operator ordering here corrects the form in (67) of [41]. The resulting stochastic description is 
unchanged. 
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This follows from operator identities of the form 

{G<-'(u,v,Lc)0(u-v/2)}ps =Ps{G<+'(u,v,Lc)0(u-v/2)}, (50) 

{G<+'(u,v,-ic)0t(u-v/2)}p,-p,{G<-'(u,v,-Lc)0t(u-v/2)}, (51) 

{M(u,v,-Lc)t/(u-v/2)}ps =Ps{M(u,v,Lc)C/(u-v/2)}, (52) 

which are derived using ( 145 > . ( I46t . ( I49> and the commutators [(/)(x), iVc] = '/'(x), and 
[C/(x),iVc] =0. 

3.6. Approximate master equation 

In the high temperature regime the master equation can be linearised to obtain a master 
equation which can, within approximations that we discuss below, be mapped to a Fokker- 
Planck equation with positive definite diffusion matrix. 



3.6.1. Approximate treatment of growth Since the dependence of the amplitude ( I38t on 
V is very sharply peaked around v = 0, we may treat the growth as a local process and 
approximate the field operators by (f){vL ± v/2) k, 0(u). We make the further approximation 
that 

G'+>(u,v,e)«G<+'(u,v,0). (53) 

This represents the fundamental result of the Quantum Kinetic Theory of Bose-Einstein 
condensation: stimulated growth of the condensate involves collisions between two 
noncondensate band atoms at right angles such that one of the atoms takes all of the energy 
of the colliding pair. The remaining atom passes into the condensate. This means that the 
growth amplitude (I38> can be approximated by its value at zero energy. From the point of 
view of [30] this amounts to neglecting the local condensate energy relative to the energy of 
the collision. 

Linearizing the forward-backward relation (I45> then gives 

G(-'(u,v,e)« (^l-^^G(+'(u,v,0) (54) 

This requires that we are in the high temperature regime in the sense that {ji — Lc)4>{'^) /ksT 
is always small. In other words the difference between the chemical potentials of the two 
bands must be small relative to the temperature. 

Combining the two approximations (I53> and ( I54t . we define the growth amplitude 

G(x) = y"d3vG(+'(u,v,0). (55) 

The effect of rotation on the growth is merely to dilate the thermal cloud, so that we can obtain 
the approximate high temperature form by using the approach of [52] to get 

G(x) = exp [2(/i - Veff (x))/fcBr] (56) 

where the rotating frame effective potential is defined in (I32> . The result (I56> follows 
immediately from the high temperature calculation of in [52] since the derivation used 
there is explicitly local, so that the chemical potential of the bath can be replaced by the local 
effective chemical potential of the rotating cloud. In fact, a more accurate form for the growth 
amplitude that includes the effects of quantum statistics and the cutoff on the noncondensate 
band distribution, found for the non-rotating case by Davis [54], can be obtained for the 
rotating scenario in exactly the same way. Although it is important to include these effects for 
certain scenarios, in this paper we restrict our attention to the approximate form i56\ . 
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3.6.2. Approximate treatment of scattering Expanding the forward-backward relation J46t 
to lowest order in e/kgT gives 

M(u,v,e)« (^l__i_^M(u,v,0). (57) 

If the condensate band is always reasonably close to equilibrium with the thermal cloud so that 
the growth linearisation is vaUd, the approximation for the scattering is the most restrictive. 
In terms of the eigenoperators of He defined by [52] 

[i/c,^m(x)] =-e™X„(x), (58) 

we have 

Lc^m(x) = e™X™(x), (59) 

and the condition ( I57> is a requirement of high temperature in the sense that 

|e„| « 2fcsT, (60) 

for the eigenvalues of Lc associated with the eigenoperators that contribute significantly to 
the expansion 

0(x)=^X™(x). (61) 

m 

The highest value of e„j available is of the order « En, so we require 

Er « 2kBT. (62) 

In [55] it was shown that the single particle energy levels of the noncondensate band are very 
closely approximated by those of the trapping potential, provided the cut-off satisfies 

2.5^i < Er. (63) 
Under these conditions, we require 

^Ji < keT (64) 
for the linearisation of the scattering to be valid. 

3.6.3. Scattering amplitude We can compute the amplitude by considering the Fourier 
transform 



M(x,k,0) = ^ j d\ e-*-^M(x,v,0) 



27/2 r f 

= Tr^ '^'^i / rf'K2F(x,Ki)[l + ^^(x,K2)] 

X 5(Ki - K2 - k)<5(co(x, Ki) - t^(x, K2)), (65) 

with a;(x, K) given by ( l3ll . 

After evaluating the momentum delta function this becomes 

,2 



M(x, k, 0) = —5-3 / i^(x, K) [1 + ^^(x, K)] 

|^2|k- X x| -k-k^) . (66) 



INC 

X 6 



Using the notation fta)(x, K) = eK-mf2xx/;i(x), where 

eK(x) = ^— + yeff(x), (67) 
2m 
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we can change to the integration variable q = K — x x) /h, and obtain 

M(x,k,0) = — / d^q frM^m + /T,p(eq(x))] 

X 5 (2q • k - k^) , (68) 

where 

= e(^-^)/LT_i - ^69) 

Note that ( I68t differs from the corresponding expression (88) of [41] only in the replacement 
V{x.) Vcb{x.), and the change of integration region to eq(x) > Ej^. A straightforward 
calculation along the lines given in [41] leads to the form 

where the high temperature condition and the fact that the operators upon which M(x, k, 0) 
acts are restricted to the condensate band have been used to obtain the final form. Note that 
the ease with which the laboratory frame calculation is shifted to the rotating frame has its 
origin in the fact that both the distributions and the edge of the integration region only depend 
on energy. In the semi-classical description we are using the energies in the two frames only 
differ by a position dependent shift in momentum and a change of the effective potential. 
In the spatial variables we find 

and the action of the scattering operator is 

d^x' M (^^,x-x',o) /(x') = ^/(x). (72) 
It is clear that the scattering term is entirely non-local. 

3.6.4. The Master equation Using the above approximations leads to the Master equation 

PC=PC\ Ha,„ + PC I g,„^th + PC I scatt (^3) 



where the Hamiltonian evolution is given by ( I23> and the growth and scattering terms are 



PC 



growth 



d^x G(x){ [[0(x),pc],'/'^(x)] + [0(x),[pc,0^(x)]] } 

d'^J^{W- Lc)cp{^)} PcAH^)] 

+ [0(x),pc{(M + ic)0^(x)}] } (74) 

and 

PcL,,tt= " Z"^'"/ ^'vAf(u,v,0){[C/(u + v/2),[C/(u-v/2),pc]] 

+ ^ [C/(u + v/2), [{LcU{n - v/2)), pcU }, (75) 

and where [ , ]+ denotes the anti-commutator The scattering terms are a non-local form of 
the quantum Brownian motion master equation [45]. 
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3.6.5. Non-Lindblad terms The Lindblad property is a highly desirable property for 
Markovian master equations since it guarantees that the positive semi-definiteness of density 
operators is preserved by the dissipation. The scattering terms in the master equation, given 
by il5\ . cannot be cast in Lindblad form. There are also growth terms which are not of 
Lindblad form, arising from the kinetic and interaction terms in the second line of ilAl . The 
kinetic terms take the form 

PC = [ d'x ^ ^ [2V(/.(x)pc (x) - V0t (x) V</)(x)pc - Pc V0t (x) V0(x)] 

kin J KbI 2m '■ ■' 

+ / '^'^ 2^;^^^W • ['/'WpcV</>t(x) + V</.(x)pc0nx) 

- pcV</)t(x)</)(x) - 0t(x)v</.(x)pc] . (76) 

If the thermal cloud varies slowly over the condensate region we may neglect VG(x) and 
recover a dissipative kinetic term in Lindblad form given by the first line of i76\ . 
The interaction terms in (I74> are 

PC = /d3x^«([0t(x)0(x)0(x)pc,0t(x)] + [0(x),pc0t(x)0t(x)0(x)]), (77) 

int J Kb i V / 

which is not of Lindblad form. However, Munro and Gardiner [56] have shown that the effect 
of non-Lindblad dissipation is usually to generate small initial transients in the time evolution, 
so we may be confident that the dissipative evolution will still give physically meaningful 
predictions. 

4. The Wigner phase-space representation 

4.1. Wigner distribution 

We expand the projected field operator in the single particle basis 

(/.(x) =^a„r„(x) (78) 

n 

where the bar notation indicates that the summation is carried out up to the energy cut-off of 
the condensate band, and 



[an,al^\ = Snm, [al,al^\ ^ [al,a'l^\ ^0. (79) 
Following [45], we may define the symmetrically ordered quantum characteristic function 



Xwi{K,K}) = jpc exp Amflt^ - A,*„a„^ |. 



(80) 



The multimode Wigner function is then given by 

j2i 



W({a„,<})= /'[]l^exp(^A;a,-A,a;]xw({A„A;}). (81) 

Moments of the Wigner distribution give symmetrically ordered operator averages, for 
example 

/T-r on / a.la„ + a„at, \ 

n d^aM'W{{a^, aU) = ( ^ % ' ^ ) • (82) 
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4.2. Operator correspondences 

Defining the projected stochastic field q;(x) and the projected derivative operators as 

a(x) EE^a„y„(x), (83) 

11 



(5a (x) " dun ' 

one then finds functional operator correspondences for the Wigner function [45] 

0(x)p ^ fa(x) + i^A^W, (86) 



2fo*(x) 

^t(xV^(.*(x)-i^)pF, (87) 

"W-^J^)^' '''' 

p<f>'^ (x) ^ (^a* (x) + ^ (89) 

which are used to map the master equation (I73> to an equation of motion for W . 
The field average corresponding to J82t is 

n / rf^a„ |a(x)Piy({a.,<}) = ^ ^nx)0(x) + 0(x)0t(x) 

= (0t(x)^(x)) + M|li^, (90) 

where the condensate band delta function 

^ X^y„(x)r„*(x') = [0(x),(/)t(x')] (91) 

n 

is always well defined at x = x' for a finite cut-off energy. The delta-function contribution 
in ( I90> exhibits a well known issue with using the Wigner representation of quantum field 
theory. The theory is ultraviolet divergent, a feature which is absent from the positive-P 
formulation [36]. Nevertheless, the Wigner representation has some utility when treating 
ultra-cold Bose gases subject to the cold-collision interaction Hamiltonian (0). Within the 
Wigner representation there is an approximation available which is known as the classical 
field method or truncated Wigner method. 

4.3. Truncated Wigner method 

By way of the operator correspondences (I86> - ( I89> . the Hamiltonian terms ilQ\ generate the 
time evolution 
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2y72 



ihVl • (x X V) + y(x) + w[|a(x)|2 - fc(x,x)])a*(x) 



?i 5a* (x) V 2m 



4n5a(x)fo*(x) V"*^"^J^ " "("^M^j ^''^ 

The third order terms present a significant difficulty in the sense that there is no equivalent 
diffusion process which would admit a formulation in terms of stochastic differential 
equations [45,57]. Although stochastic difference equations can be found which are formally 
equivalent to the generalized Fokker-Planck equation [58], they are difficult to use. The 
truncated Wigner method [36,44,59-62], is based on the observation that the third order 
derivatives are small when the modes of interest are highly occupied throughout. If we neglect 
the third order terms in i92\ . the resulting equation of motion reduces to a Liouville equation 
since it lacks any diffusion terms, and the equation of motion for the classical field q;(x) 
has the superficial appearance of the Gross-Pitaevskii equation. It should be stressed that 
there are two important differences, (i) As we have seen, moments of the Wigner function 
give symmetrically ordered operator averages which must be calculated over a large number 
of trajectories, (ii) The initial state must also be chosen according to the Wigner transform 
of the initial density matrix. Although the equation takes the form of the Gross-Pitaevskii 
equation, only in the limit that the initial Wigner function is a Dirac delta distribution do we 
recover the Gross-Pitaevskii theory. We note Polkovnikov has shown that truncated Wigner 
approximation arises naturally as the first in a hierarchy of quantum corrections to the classical 
field evolution described by the Gross-Pitaevskii equation [63, 64]. 

The application of the truncated Wigner approximation to the generalized Fokker-Planck 
equation arising from ( 17 3> leads to a genuine Fokker-Planck equation with a positive definite 
diffusion matrix which can then be mapped to a stochastic differential equation of motion for 
the condensate band. 

5. Condensate band Fokker-Planck equation 

We can now use standard methods [45] to map the Master equation to a genuine Fokker-Planck 
equation for the Wigner distribution. 

5.7. Hamiltonian tenuis 

Defining the Gross-Pitaevskii Hamiltonian corresponding to ilO\ 

iJcp = / a*(x) f — — + iMl • (x X V) + y(x) j a(x) 

+ ^ J d^x a*(x)a*(x)a(x)a(x), (93) 

the GP operator Lqp 

VLopaix) EE (94) 
and neglecting third order derivatives in ( I92> . we can express the Hamiltonian term as 



dW 
~dt 



d^x \ ( -^LGpa(x) ) + c.c. \ W. (95) 

Ham J [ da[x) \ n 
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Note that Lqp is implicitly projected since the action of a projected functional derivative is 
equivalent to non-projected functional differentiation of a projected operator - a consequence 
of the Hermiticity of projectors. We can express \95\ as 

dW 
~dt 

from which can use the usual correspondence rules of functional Fokker-Planck equations [36] 
to identify the drift term as the projected Gross-Pitaevskii equation of motion [65] for a(x, t), 
as seen in the first line of ( I106> . 



= I d"^ { \^--rLGpa{^) ] + c.c. } W, (96) 



5.2. Growth terms 
The master equation terms ilAl give 
dW 



J d^xG(x)< 



in - LGp)a(x) 



growth J I Sa[x) \ ksT 

6 S 1 

Ja*(x) (5a(x) J ^ ^ 

which contributes a drift term and a thermal noise term to the stochastic Gross-Pitaevskii 
equation. 

5.3. Scattering terms 

Our development of the formalism thus far has largely followed [41], the only modification 
being the — • L term in the single particle Hamiltonian arising from the rotation. Similarly, 
for the scattering there are extra terms arising from the formulation in a rotating reference 
frame. 

To derive the scattering terms in the Fokker-Planck equation we need to evaluate (175 > . It 
is easily seen that 

LcC/(x) = -fW-jc(x) (98) 
where the condensate band current operator 

jc(x) = ^ ([V0t(x)]0(x) - </)t(x)V0(x)) -{nx x)(/.t(x)<^(x) (99) 

exhibits a centrifugal term in the rotating frame. The corresponding Gross-Pitaevskii current 

jGp(x) = -^[a(x)Va*(x) - a*(x)Va(x)] - (Ox x)a*(x)a(x) (100) 

arises in the drift vector of the Fokker-Planck equation, which takes the form 
dW 
'dT 



d^x I d-^x' M I — - — , x' - x, 



S , iW • jGp(x) 
— a(x ) + c.c. 




(101) 
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5.4. Full Fokker-Planck equation and its stationary solution 

The Wigner distribution evolves according to the Fokker-Planck equation 



dW _ dW 
~dt ~~dt 



dW 
~dt 



growth 



dW 
~dt 



By defining the GP number operator corresponding to ( I47t 

iVcp = j d'^x a* (x)a(x), 
the equation of motion can be written as 



dW 
'df 



d'^x 



iW SHgp 



Sa(x) [ h (5q!*(x) 
S r iW SHgp 



Sa*{x) [ h (5a(x) 



G(x 
-G(x 



SW 
Sa* (x) 

SW 



W 5{iiNgp-Hgp) 



ksT 5a* (x) 
W S{fiNGP-HGP 



d'x I (fx M 

SW 



X <^a*(x') 



W SHgp 



5a (x) 
a(x) 

-a(x') 



5 



(5a(x) 



5a(x) 
SW 



a*(x) 



(5a*(x)_ 

W SHgp 

5a(x') ksT Sa{x') 
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(102) 



(103) 



.(104) 



_(5a*(x') fci3rfo*(x') 

Irrespective of the forms of G(x) and M(u, v, 0), the Fokker-Planck equation has the grand 
canonical stationary solution 



Ti/ '^j.Ngp-Hgp\ 
l^ocexpl ■ 



(105) 



6. Rotating stochastic Gross-Pitaevskii equation 

6.1. Full form of the SGPE 

In the rotating frame the full non-local form is given by the stochastic differential equation in 
Stratonovich form 

(S)da(x,t) = - ^VLGpa{yi)dt 

h 



G(x) 
keT 



{p - LGp)a{x.)dt + dWa (x, t) 



d^x' M 



■,x-x',0) -^V-jGp(x')rfi 



-I- idWMi^, t)a(x) 



(106) 



where the growth noise is complex, the scattering noise is real, they are independent of each 
other, and satisfy 

dWa{-^,t)dWG{-x! ,t) = 2G(x)(5c(x,x')dt, (107) 
dWG{^,t)dWG{y^ ,t) dVF5(x,t)dW^5(x',i) = (108) 



dWui^, t)dWM{^, t) = 2M ( X - x', ) dt. 



(109) 
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6.2. Approximate non-local form 

In the rotating frame the approximate non-local form of the SGPE [41], which is probably the 
most practical form for realistic simulations, is given by the Stratonovich SDE 

(5')da(x,t) = - lrLGpa{x.)dt 

h 

+ -1^^ ■ JGP (x)di + idWM (x, i)a(x) I , (110) 

where the growth noise is complex, the scattering noise is real, they are independent of each 
other, and satisfy 

dW^{-K,t)dWG{yi\t) = 2G(x)5c(x,x')dt, (111) 
dWG{^,t)dWG{-x! ,t) = dVK<5(x,i)dVF5(x',t) = (112) 

dWM{^.t)dWM{:^ ,t) = -^^5c(x,x')di. (113) 

V — 

There are some corrections to mention here in the scattering terms of the forms ( I106> and 
Jl 10> . compared to equations (118) and (122) of [41]. In particular, we note that the scattering 
term in [41], proportional to a*(x, t)Lca(x, t) — Q!(x, t)(ic'Q!(x, t))*, is subtly different to 
the corresponding expression in jl06> . M 10> . V • Jgp, since there are no projectors in the latter. 
Furthermore, the form of M arising in our ( I70> and (^Oj is bigger than the form defined in 
equation (97) of [41] by a factor of (27r)3. 

7. Theory of irreversible vortex nucleation 

Now that we have developed a theory that includes the effects of a rotating thermal cloud 
on a trapped BEC, we will examine the connection between a simplified form of the theory 
described in section|6land the other approaches that have been used to investigate the dynamics 
of vortex lattice formation. 

Much effort has been directed at understanding the process of vortex lattice formation 
in trapped Bose-Einstein condensates [29,33,66,67]. It is well known that the minimum 
energy state of a BEC in a rotating reference frame usually contains vorticesf . An obvious 
consequence of this fact is that there must be some irreversible dynamics for the system to 
proceed from a laboratory frame ground state to a rotating frame ground state; physically this 
evolution arises from either (i) contact of the condensate with a rotating thermal cloud, or (ii) 
turbulent motion arising from dynamical instability. Mechanism (ii) has been already been 
the subject of significant research [1 1, 33, 34, 68]. In this paper we are concerned with (i). 

Although much work has been carried out using phenomenological and hydrodynamic 
theories [29,66,67], the precise description of the thermodynamic scenario has, until now, 
been an oustanding unsolved problem. Furthermore, it is to be expected that the combination 
of dynamical instability and nonlinear interactions generates a thermal component, which 
in turn acts as a rotating reservoir that imparts angular momentum to the condensate. We 
therefore wish to formulate a precise description of a trapped condensate in contact with a 
rotating thermal cloud, since it is this aspect of the theory that provides the irreversibility 

t In the regime of slow rotation tliere are states tliat have nontrivial phase variation but are nevertheless vortex 
free [50]. 
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required for vortex lattice formation. As is often the case in the theory of dilute gases, we are 
able to develop a first principles formulation of the theory that goes beyond phenomenology. 

The recent history of the phenomenology of damped condensates can be traced to the 
work of Choi, Morgan and Burnett [28] who proposed that dissipation may be incorporated 
into the GPE theory by introducing a small imaginary component to the time evolution. While 
this approach has been very successful in practice, the theory is entirely phenomenological. 
The outstanding issue we have resolved is the consistent description of a rotating thermal 
reservoir, including the noise arising from Bose-enhanced scattering into the condensate band. 
We will now show that a simplified form of the SGPE theory leads in a natural way to a picture 
of GPE evolution in 'complex time', as is to be expected on physical grounds for a system 
coupled to a thermal reservoir The damped GPE approach of [28,29,41,66-68] is very close 
to this kind of description; here we make the connection rigorous. 



7.1. Vortex growth equation 

An equation of motion which describes the dissipative dynamics of vortex lattice formation 
can be obtained from ( II 10> with some additional approximations. The scattering terms are 
neglected on the grounds that they primarily lead to a renormalisation of the effective growth 
rate [24,55], and so can usually be regarded as a second order effect. We make the simplifying 
approximation that the thermal cloud is spatially smooth over the region of the condensate, 
and then set dWc ^ 0, a ^ hG{x)/kBT 7, and P ^ 1, to find 



= icpV-lx, t) + i7(m - Lgp)V'(x, t). (114) 

ot 

This is the simplest description of a condensate in contact with a noiseless rotating thermal 
cloud at fixed chemical potential /i. We make the further transformation — ae^*^*/'*, to 
obtain 

■^da^^^^_^^^^^^_ (115) 
ot 

and it is clear that the condensate wavefunction is evolving in 'complex time' obtained from 
the standard GPE evolution by the replacement dt dt{l — ij). From ( I93> and ( I103> it is 
easy to show that 

Q-^ = — ^ / d-'x |(/i-LGp)a(x,t)| , (116) 

which shows that the evolution minimises Hqp — iiNqp and the wavefunction evolves into 
the rotating frame ground state of the GPE with the same chemical potential as the thermal 
cloud - a vortex latticef . 

Starting from first principles we have arrived at a satisfactory theory of the dissipative 
dynamics of vortex lattice formation due to energetic (thermodynamic) instability - without 
thermal noise. In essence the theory is obtained by shifting to the cloud frame, and then 
carefully accounting for interactions between cloud and condensate. The principal aim of this 
paper is to provide a theory which also accounts for the physical role of thermal noise in this 
process. We now briefly outline the differences between this and previous formulations. 

f For simplicity we have set 'P — > 1. However, note that <1 161 still holds if we relax this condition [69]. In other 
words the evolution toward the ground state expressed by <1 161 also occurs when the theory explicitly includes a high 
energy cutoff separating condensate and thermal cloud. 
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7.2. Phenomenological and hydrodynamic models of dissipation 

The numerical vortex nucleation studies of Tsubota et al [29], and later Kasamatsu et al [68], 
were based upon a phenomenological damped Gross-Pitaevskii equation. This description 
was introduced in [29] on the grounds that it gives the correct equilibrium solution, and that 
it follows naturally from the damped GPE approach of [28]. 

The physical argument for such an equation, was first provided by Gardiner et al [30], 
and the equation of motion, found via reasoning we will explain below, was subsequently used 
by Penckwitt et al [66] to investigate vortex lattice formation. Like the Tsubota equation, 
the derivation involved a hydrodynamic approximation, in particular it made use of the local 
energy concept and the approximation /ic(x, i)^' ~ ihdip/dt, where /ic(x, i) denotes the 
local condensate chemical potential. As described later in [68], the equation of motion may 
also be obtained via similar reasoning from the generalized finite temperature GPE derived 
by Zaremba et al (ZNG) [70]; we briefly reiterate the argument here. 

In the rotating frame the ZNG equation takes the form 

_dip{yi,t) _ ( h^V^ 



'^^t " + ■^trap(x) - nL, + «|i^(x,t)|^ + 2un(x,i) - iTj yj{x,t), (117) 

where n(x, t) is the noncondensate density, and F describes collisions with the 
noncondensate. This equation of motion has been derived within a hydrodynamic 
approximation, utilising the notion of local energy attributed to a smoothly varying condensate 
wavefunction. We further neglect the noncondensate density under the assumption ^ h, 
and note that under a local equilibrium distribution for the thermal cloud the growth rate is 
proportional to the difference between the local chemical potentials of the condensate and 
noncondensate F cx /inc(x, t) — /Xc(x, t); the constant of proportionality will be denoted by 

7- 

We now make the crucial approximation that generates the dissipative GPE evolution, by 
putting /Xc(x, ~ —ih{dip/dt). To reach the final form, one then proceeds by neglecting 
the space and time dependence of the noncondensate chemical potential (fiad^, t) A'nc — 
jj), and making the further transformation — > ?/)e~*'^*/''. We finally arrive at the equation of 
motion introduced by Tsubota et al [29] 

{^ ~ = ("^ + ^*-pW ~ + t)P - m) ^(x, t). (118) 

Casting this in the form 

^^^^^ 

ot 1 + 7 

we see that in the small 7 limit we recover the vortex growth equation ( I115> . However, for 
non-negligible 7 the imaginary and real time dynamics generated by ( I118t are slower than 
those of the vortex growth equation il 14> . a consequence of shifting to complex time without 
accounting for the change in the length of the time vector in the complex plane. 

We stress that neither of these approximations have been made in obtaining the SGPE 
([^roj, and the simpler vortex growth equation il 14t . Furthermore, the simple growth equation 
has all the necessary physical properties for the correct damped GPE description of vortex 
lattice formation. 



8. Conclusion and Outlook 



We have extended the SGPE theory developed in references [30,41] to include the possibility 
that the thermal cloud may be in a state of rotation - a situation of considerable experimental 
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interest. The theory is suitable for treating the experimentally interesting regime of rotating 
Bose-Einstein condensation [27], where one may expect that vortices can form during 
condensation, effectively being 'pinned' in the condensate so that subsequent scattering into 
vortex states becomes Bose-enhanced. This is an entirely different physical regime to the 
coherent vortex formation mechanism arising from stirring a pre-existing condensate [9]. 

We have examined the connection between the SGPE theory derived here and the 
phenomenological and hydrodynamic approaches that have been used thus far to model vortex 
lattice formation. We have shown that a simplified vortex growth equation has all the desirable 
properties of the other mean field approaches, while having a wider range of validity and 
requiring much less severe approximations in its treatment of the condensate evolution. 

The theory presented here contains two additional features that have not yet been 
investigated in the literature on vortex lattice formation in Bose-Einstein condensates: thermal 
noise and number conserving scattering. While we have not discussed the influence of 
scattering on vortex dynamics here, one may expect that vortices will act as strong scatterers, 
and indeed that this process may be far more significant for vortex lattices than the weak 
renormalisation of the growth rate observed for vortex free condensates [24]. Thermal noise 
is an essential feature of any dissipative theory, and its effect on vortex lattice formation, both 
for the mechanical stirring and defect pinning nucleation mechanisms, will be the subject of 
future work. 
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